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ABSTRACT 

We compare the prf/and its low-order moments (variance and skewness) of the 
smoothed PSCz galaxy density field and of the corresponding simulated PSCz 
look-alikes, generated from N-body simulations of 6 different DM models; four 
structure normalized with F — 0.25 and cg = 0.55il^°-^, one COBE normal- 
ized and the old SCDM. The galaxy distributions are smoothed with a Gaus- 
sian window at three different smoothing scales Rsm = 5, 10 and 15 ft,~^Mpc. 
We find that the simulation PSCz look-alike pdfs are sensitive only on the 
normalization of the power spectrum, probably due to shape similarity of the 
simulated galaxy power spectrum on the relevant scales. We find that the only 
models that are consistent, at a high significance level, with the observed PSCz 
pdfaie models with a relatively low power spectrum normalization (erg = 0.83). 
From the phenomenologically derived erg ~ moments relation, fitted from the 
simulation data, we find that the PSCz moments suggest cg ~ 0.7. 

Keywords: cosmology: theory - galaxies: general - large-scale structure 
of universe - Infrared: galaxies 



1 INTRODUCTION 

The study of the galaxy or cluster one-point probability 
density function (pdf) , defined as the probability of hav- 
ing A'' cells of volume V with overdensity S{= p/ (p) — 1), 
has been used in order to probe the Gaussian or not na- 
ture of the cosmic density field as well as cosmological 
probe, since the pdf depends on the two-point correla- 
tion function or its Fourier Transform, the power spec- 
trum, but also on higher order correlations. 

The pdf oi the primordial density field evolves due 
to gravitational instability and even if it follows Gaus- 
sian statistics originally, it eventually becomes non- 
Gaussian by the mere effects of gravity. One therefore 
may envision that smoothing on large enough scales the 
present day galaxy density field one can recover the ini- 
tial conditions. This procedure is unstable due to the 
fact (1) that galaxies are probably a biased tracer of 
the underline density field and (2) the central limit theo- 
rem ensures that heavy smoothing will Gaussianize even 
an inherently non- Gaussian field. However, elaborate in- 
version techniques have been devised and applied with 
varying level of success to real and simulation data (see 
Monaco et al. 2000 and references therein) and the gen- 
eral outcome is that indeed the present-day galaxy dis- 
tribution appears to originate from the effects of grav- 
itational instability on an initial Gaussian fiuctuation 
field. 

Many authors studying the pdf statistics of differ- 



ent samples of extragalactic objects, assuming linear 
biasing, have claimed that if the large-scale distribu- 
tion of galaxies trace the underlying matter distribu- 
tion then the observed pdf can be approximated us- 
ing Gaussian initial conditions (cf. Bouchet et al. 1993; 
Gaztanaga & Yokohama 1993; Sheth et al. 1994; Kof- 
man et al. 1994). Other authors, using N-body simula- 
tions, studied the evolution of the pdf and its moments 
(cf. Coles & Frenk 1991; Juszkiewicz et al. 1993; Coles 
et al. 1993; Kofman et al. 1994; Colombi 1994) and con- 
cluded that the normalized skewness of the matter dis- 
tribution, Ss — {S'^)/{S^) depends on the shape of the 
power spectrum. 

Saunders et al. (1991) used the lower-order mo- 
ments of the pdf (the variance and skewness) of the 
sparse sampled IRAS 0.6Jy sample (QDOT) to put 
constraints on cosmological models. They found that 
on large scales they are larger than the corresponding 
moments of QDOT look-alikes, generated from the old 
standard CDM model. In this paper we will utilize the 
completed PSCz-IRAS redshift survey in a similar at- 
tempt. 

The layout of this paper is as follows: in Section 2 
we describe briefly our simulations and the considered 
models; in Section 3 we present the analysis method and 
the results obtained; we state our main conclusions in 
Section 4. 
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2 THE DATA 

2.1 The PSCz sample 

In this work we use the recently completed PSCz-IRAS 
flux-limited redshift survey which is described with de- 
tails in Saunders et al. (2000). It is beised on the IRAS 
Point Source Catalog and contains ~ 15500 galaxies 
with flux Slim > 0.6 Jy covering the 84% of the sky. 

Due to the fact that the PSCz catalog is a flux- 
limited sample there is the well-known degradation of 
sampling as a function of distance from the observer 
which is represented by the selection function. It is nec- 
essary therefore to recover the true galaxy density field 
from the observed flux-limited sample. This is done by 
weighting each galaxy by <p^^{r), where the selection 
function, (^(r), is defined as the fraction of the galaxy 
number density that is observed above the flux limit at 
some distance r In this work we use a luminosity func- 
tion of the form assumed by Saunders et al. (1990) with 
L* = 10®-*^ h^LQ, G = 0.711, a = 1.09 and C = 0.0308 
(cf. Rowan- Robinson et al. 2000). 

In order to take into account (1) the 16% of the sky 
which is devoided of data, due to high cirrus emission 
areas and unobserved regions, and (2) the distortion of 
galaxy positions due to peculiar velocities, we utilize 
the corrected galaxy catalog of Branchini et al. (1999). 
These authors followed the Yahil et al. (1991) prescrip- 
tion to flU the galactic plane with synthetic objects that 
reproduce the mean density of galaxies of the nearby 
regions while at higher galactic latitudes they fill ran- 
domly the masked regions with synthetic objects, repro- 
ducing the overall mean number density of PSCz galcix- 
ies. Furthermore, in order for them to minimize the so 
called "Kaiser" effect (Kaiser 1987) and recover the true 
3D galaxy distances they used an iterative technique 
and reconstructed the PSCz galaxy density field for dif- 
ferent values of /3 = f2o'®/6iRAS. In our present analy- 
sis we use their (3 = 0.5 reconstruction, since recently 
Branchini et al. (2001), using the VELMOD analysis 
to correlate the PSCz density field with the SFI galaxy 
velocities (cf. Haynes et al. 1999) found Piras — 0.42. 



2.2 N-body Simulations 

We use mock PSCz catalogs, generated from six large 
cosmological N-body simulations (see Cole et al. 1998), 
in order to investigate whether the pdf and its moments 
can discriminate between models. The PSCz look-alikes 
have been generated by E. Branchini, choosing suitable 
LG-like observers (having similar to the observed Local 
Group velocity, shear and overdensity) , and he gener- 
ously provided them to us. Note that due to our igno- 
rance in assigning galaxy formation sites to the DM halo 
distribution and for a consistent treatment of all mod- 
els, a biasing factor of 1 has been used in generating the 
PSCz look-alikes. 

Four of the models that we consider are structure 
normalized with F = 0.25 and fluctuation amplitude in 
8 /i-^Mpc scale of cts = 0.550"" '^ (Eke, Cole, & Frenk 
1996), which are: (1) a flat low-density CDM model 
with Q,o = 1 — Ha = 0.3 (Aqdmi) (2) a model with 



Table 1. probability of significant differences between 
the model-model or model-data pair pdfs. Each probability 
value corresponds to the pair formed between the indicated 
model or data in the first column and the first row. In bold we 

indicate the results of the PSCz-model pdf comparison. For 
example, a high probability value (say > 0.1) means that 
the corresponding two pdfs show no significant differences 
between them. When the probability is less than 10~^ we 
write 0.00. 
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57o = 1— f^A = 0.5 (AcDM2) (3) an Slo = 1 model (tcdm) 
and (4) an open model with Qo = 0.5 (Ocdm); a COBE 
normalized Einstein de Sitter model with erg = 1.35 
(Ccdm) and flnally the old 'standard' CDM model with 
r = 0.5 and ag = 0.55 (Scdm). 

For each cosmological model we average results over 
10 nearly independent mock PSCz catalogs which are 
treated in exactly the same way as the real PSCz data. 



3 THE PDF DERIVATION 

We obtain a contirmous galaxy density held by smooth- 
ing the discrete distribution of PSCz galaxies on a 
grid. We use a Gaussian kernel with smoothing radius. 
Rem, varying from 5 to 15 Mpc. The grid size that 
we use is set equal to Rsm- If p(x) is the smoothed den- 
sity field at a grid point, then the relative fluctuations 
are given by (5(x) = p(x)/(p) — 1, where the average 
density (p) does not depend on Rsm- The frequency 
distribution of S's is the pdf. 

Basilakos, Plionis & Rowan-Robinson (2001) in- 
vestigated the systematic biases that enter in such a 
smoothing procedure due to the convolution between 
the PSCz selection function and the constant smooth- 
ing radius. The effect is such that the derived smooth 
density held, at large distances, is overestimated in re- 
gions where galaxies are detected (due to the large 0(r) 
weighting), while it is underestimated in underdense re- 
gions. They devised a phenomenological procedure that 
corrects effectively the high density regions for such ef- 
fects. Since this procedure is multiplicative in nature, it 
is not as effective in recovering the underdense regions 
and thus we confine our analysis within a limited radius. 



PSCz pdf 3 




rmax (= 80, 110 and 120 h'^ Mpc for the R,m = 5, 10 
and 15 Mpc fields, respectively), within which we 
have verified that the derived density field is mostly un- 
alTected by the above mentioned bias. This was achieved 
by comparing the pdf and its moments of PSCz look- 
alikes and the full 3D "galaxy" distribution in the sim- 
ulations. Even within rmax the raw values of the vari- 
ance and skewness are overestimated by ~ 20%, but 
only for the Rsm = 10 and 15 Mpc fields. Applying 
the Basilakos et al. (2001) phenomenological correction 
procedure does recover the correct values of these mo- 
ments. 

In order to investigate how discriminative between 
models is the pdf, we have performed a statistical 
test, comparing the model pdf between them. We find 
that the present pdf analysis can discriminate only dif- 
ferent (Tg models, irrespectively of the DM content and 
the geometry. By using the whole 3D simulation data, 
we have verified that neither shot-noise nor the PSCz se- 
lection function is responsible for such behaviour; which 
is probably due to the shape similarity of the derived 
"galaxy" power spectrum on the relevant scales (see fig- 
ure 2 of Cole et al. 1998). This is particularly true for 
the OcDM and Acdm2 models that both have F = 0.25 
and erg — 0.83. 

Furthermore, we find that the most effective, in 
discriminating power, smoothing scale is the Rsm ~ 5 
h-^ Mpc one. In Table 1 we present the probabilities, 



Pj,^2 , of significant differences between any two pdfs for 
the two indicated smoothing scales. It is clear that the 
statistical test identifies 3 families of consistent among 
them models, all of which are based on their relevant erg 
normalization. 

In figure 1 we present a comparison between the 
PSCz pdf and those of the 3 families of models. From 
this figure as well as from table 1, it is clear that the 
models with erg — 0.83 are the only ones consistent with 
the data. 

3.1 Moments of the pdf 

The variance and the skewness 7 are defined as the 
second- and third- order moments of the S field, respec- 
tively. We therefore have 

n 

71 

7=i5]5'(x), (2) 

where n is the total number of grid points within rmax ■ 
Gaztanaga & Yokoyama (1993) have shown that the 
smoothing process itself suppresses considerably the 
shot-noise effects. A proper treatment of these effects 
should entail subtracting the "Poisson" terms from a 
and 7 but in conjunction with the effects of smooth- 
ing, which is not a trivial task. However, since all the 
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Table 2. Moments of the smoothed density distribution of 
the mock and real PSCz data. The uncertainties of the mod- 
els are estimated from 10 PSCz look-alike realizations of each 
model while those of the PSCz data are estimated by dividing 
the available volume into 4 independent regions. 
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mock PSCz distributions have the same (p) and we treat 
them similarly, the possible shot-noise effects are rela- 
tively and qualitatively cancelled out in the model and 
data intercomparison. 

In Figure 2 we plot the plane, following Coles 
& Frenk (1991), for the PSCz density field (red filled 
dots) and for the simulation mock PSCz fields. The re- 
sults of the Rsm ~ 5, 10 and 15 Mpc field are 
represented as open circles, triangles and squares, re- 
spectively. Each simulation point represents one of the 
10 realizations, so that the scatter represents the effect 
of cosmic variance. The PSCz errorbars represent the 
scatter estimated by dividing the whole available vol- 
ume in 4 equal and independent subsamples. 

We can estimate the reduced skewness by least- 
square fitting a relation: log 7 = log((T^)° +log S3, which 
for a = 2 is a generic prediction of perturbation the- 
ory, of the hierarchical and log-normal clustering mod- 
els on linear scales (cf. Coles & Frenk 1991). According 
to Bouchet et al. (1992) the reduced moments are in- 
dependent of the density and geometry of cosmological 
models and for this reason they have been used to dis- 
criminate between different biasing models rather than 
cosmological models (cf. Gaztanaga & Frieman 1994). 
Indeed, keeping a = 2 we find S3 ~ 2.7 for all the mod- 
els, with scatter over the 10 realizations, of SS^ ~ 0.4. 
This value is in agreement, within the errors, with the 
PSCz value (Sf ■^'^^ ~ 2.1 ± 0.6), where the scatter is 
from the 4 independent volumes used in the estimate. 
Furthermore, this value is in excellent agreement with 
S3 ~ 2.2, estimated from the IRAS 1.2Jy sample (cf. Fry 
& Gaztanaga 1994). Leaving unconstrained the value of 
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Figure 2. Scatter-plots for the variance-skewness relation. 
Each point refers to a single realization of the corresponding 
model. Squares, triangles and circles are for Rsm = 5, 10 
and 15/i~^Mpc respectively. Large filled symbols represent 
the real PSCz moments. 



a we find, for both PSCz data and models, which by 
construction have Gaussian initial conditions, S3 values 
in the range 3.1 — 3.6 and a ~ 2.2 ± 0.15. The deviation 
from Q = 2 is probably due to the effects of shot-noise. 

These results argue against large-scale non- 
gaussianity in the Universe since, according the study of 
Coles et al (1993), non-gaussian models show departures 
from such a, — j relation. 

In Table 2 we report the real and simulation PSCz 
pdf moments. Again from both the figure 2 and table 
2 it is evident that the PSCz results are nearer to the 
£78 — 0.83 simulations. 



3.2 Recovering the as of our Universe 

What is also evident, from figure 2, is that there is a 
tight scaling between as and the values of and 7 in 
the simulation PSCz look-alike fields. We have verified 
this by fitting a linear relation to the robust Ram = 5 
Mpc result, to find 



erg = 1.062(±0.06)(7^ 
: 0.189(±0.02)7 



0-8 



+ 0.313(±0.03) 
f 0.661(±0.03) 



with a correlation coefficient of R — 0.93. 

Since from the simulation data we have seen that 
the pd/ moments depend only on the value of as, we can 
now attempt to recover, by inserting our PSCz moments 
to the above relations, the as of our Universe. We find 

as ~ 0.7 ±0.03 , 
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where the quoted uncertainty is the scatter of the mean 
from using the raw or corrected moments (see section 
3) and the ag-variance or skewness relation. This result 
is in agreement with many recent studies (cf. Eke et 
al 1998; Einasto et al. 1999; Moscardini, Matarrese & 
Mo 2001) and in particular with the most recent clus- 
ter X-ray luminosity function evolution results from the 
ROSAT Deep Cluster Survey (Borgani et al 2001). 



4 CONCLUSIONS 

We have derived the pdf moments of the PSCz galaxy 
catalog and compared them with the corresponding mo- 
ments of mock PSCz distributions, generated from six 
different cosmological models. The pdf and its low-order 
moments were estimated from the Gaussianly smoothed 
density field on a grid, using three different smoothing 
scales, Rsm = 5, 10 and 15 /i~^Mpc. For each model we 
have used 10 PSCz-like realizations to have a handle of 
cosmic variance. We find that this type of analysis can 
discriminate models only on the basis of their different 
as values. We find that from the three erg families of 
CDM models (as = 0.55, 0.83 and > 1.1) the PSCz pdf 
and its lower order moments are consistent with models 
having intermediate values of as{= 0.83). 

Furthermore, using a phenomenological erg — 
moments relation, derived from the model PSCz look- 
alike pdfs, we find that the PSCz data strongly suggest 
erg ~ 0.7 ±0.03. 
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